Discrete light localization in one dimensional nonlinear lattices with arbitrary non 
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We model discrete spatial solitons in a periodic nonlinear medium encompassing any degree of 
transverse non locality. Making a convenient reference to a widely used material -nematic liquid 
crystals-, we derive a new form of the discrete nonlinear Schrodinger equation and find a novel family 
of discrete solitons. Such self-localized solutions in optical lattices can exist with an arbitrary degree 
of imprinted chirp and a have breathing character. We verify numerically that both local and non 
local discrete light propagation and solitons can be observed in liquid crystalline arrays. 

PACS numbers: 42.65.Jx, 42.65.Tg, 42.70.Dm 



I. INTRODUCTION 

Optical energy localization in lattices has become an 
important branch of contemporary nonlinear science, due 
to a wealth of basicphysics and potentials for light 
switching and logics 0, (^'^d references therein). 
Attention has been recently devoted to light propaga- 
tion in the frame of tunable discreteness -i.e. in lattices 
with an adjustable period, index contrast, nonlinearity. 
Examples to this extent are waveguide arrays in photore- 
fractives and droplet arrays of a Bose Einstein con- 
densate 0,0,0, 111- In nematic liquid crystals (NLC), ma- 
terials encompassing a large non resonant reorientational 
response, spectrally extended transparency, strong bire- 
fringence and mature technology , excitation as well as 
switching/steering of discrete solitons has been reported 
in voltage-tunable geometries 0, 0, 0| . Discrete soli- 
tons in NLC result from the interplay between evanescent 
coupling (owing to discreteness), molecular nonlinearity 
(leading to progressive mismatch as the extraordinary in- 
dex increases) and non locality (owing to intermolecular 
elastic forces). The latter aspect, des pite its role in sev- 
eral systems 111 111 111 111 111 111 [1II23 including NLC 
[2lll2^ . has only been discussed in the framework of ID 
discrete lattices with reference to first order contributions 
[Tsj and long range dispersive interactions 0, 0, , a 
general description of discrete solitons in the presence of 
a transverse nonlinear non locality being still lacking. 
In this paper, for the first time and with explicit ref- 
erence to a physical system of interest -i. e., nematic 
liquid crystals- we model discrete light localization in 
media with an arbitrary degree of non locality, elucidat- 
ing the interplay between non locality and nonlinearity 
on soliton dynamics. Using coupled mode theory to de- 
rive the governing equations (CMT), we demonstrate 
the existence of a new family of chirp-imprinted discrete 
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FIG. 1: Example of a discrete optical lattice in a cell of 
thickness d filled with undoped nematic liquid crystals: (a) 
front view, (b) side view. A periodically varying electric field 
Ex is applied trough an electrode-array across x and alters the 
mean molecular angular orientation 9(x, y), inducing an index 
modulation with the same period A. The top graphs sketch 
the distribution E^ versus y (left) and versus propagation z 
(right), respectively. 
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breathers which could not be sustained by a purely local 
response. Finally, we verify the theoretical predictions by 
numerical experiments with a standard NLC. The paper 
is organized as follows. Section ITU introduces a model of 
liquid crystals and carries out an original reduction to a 
non integrable discrete nonlinear non local Schrodinger 
equation, outlining the novelties with respect to previous 
studies dealing with non locality. In Sec. IIIII we apply 
a variational approach using a convenient soliton ansatz 
and derive the differential equations for the evolution of 
soliton parameters in propagation. We highlight the im- 
pact of non locality on soliton generation and demon- 
strate a novel family of discrete chirped solitary waves, 
never reported before. Finally, in Sec. IIVI we perform a 
full numerical simulation of the actual liquid crystalline 
system and demonstrate the excellent agreement with our 
analytical predictions. We conclude by emphasizing how 
the examined NLC-lattice offers the rare possibility to 
observe both local and non local light propagation in one 
and the same system. 
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II. THEORETICAL APPROACH 

We consider light propagation in a thin fihn planar 
waveguide of nematic liquid crystals, subject to a peri- 
odic transverse modulation along y and across x (Fig. 
NLC consist of rod-like molecules which, electrically po- 
larized across x, react to and reorient towards the field 
vector in order to minimize the free energy 0. Under 
"planar" anchoring conditions at top and bottom inter- 
faces of the cell, the mean angular orientation of the NLC 
molecules (i.e., the molecular director) is conveniently 
described by their angle with the axis z in the plane 
(x, z), as sketched in Fig. ^ This identifies the NLC 
optic axis with respect to the propagation wavevector 
of a light beam injected in the cell. If = ?t.| — 
is the NLC optical birefringence (with and nj_ along 
or orthogonal to the director, respectively), an electric 
field (static or low frequency) applied across x, constant 
in z and periodic along y (Fig. can reorient the di- 
rector and determine a one-dimensional optical lattice 



with index modulation n (cc, y) 



' 9{x, y) for 



e-polarized light. We assum e an applied electric field 
Ex{x,y) — Eq ■ \/^ + ^F{y), with a zero mean- value 
F{y) — F{y + A) and an arbitrary e < 1. In actual 
experiments, is determined by the bias V{x,y) oc 
x(l -I- V{y)), with V{y) — V{y -\- A) applied through an 
array of parallel finger-electrodes 

In the framework of the elastic continuum theory Q , the 
director distribution ^0(2;, y) at rest -i.e., with no injected 
light - can be obtained by minimizing the NLC energy 
functional 



{dliOoi)) « 0), we obtain: 



dy^ 9a dy dy 



3K 



2K 



(3) 



Eq. (O models the all-optical response of the NLC lat- 
tice. It can be cast in the integral form ip = J J G{(^ — 
~ y)\^{C,v)\'^d,Cdrj, with G{x,y) the Green func- 
tion. For a sufficient bias to induce an array of single- 
mode channel waveguides with nearest-neighbor cou- 
pling, CMT (in the tight-binding approximation) yields 
the z evolution of each eigenmode: 



. dQn 
dz 



' C{Qn+l+Qn-l) 

+ Qn^r,„,„|Q™|2 = (4) 



with A = J2nQn{z)fQ{x,y)exp{-i/3z), fg and /3 the 
modal eigenfunction and eigenvalue, respectively, 
the mode-power in the nth-channel, C the coupling 
strength and: 

T,n,n ^'Y' J IfQi^^y - 'n^)\'^n2,^ix,y) 

X G(C-x,r;-y)|/Q(C,r/-mA)|2dxdydCdr, (5) 

the nonlinear overlap integral, with n2,^{x,y) = 
2"(eo)''^° • Factorizing the Green function as G{x, y) = 
G{x)Gp{y)Ge{y), with Gp{y) = Gp{y + A) and the enve- 
lope Ge{y) wider than the guided mode, the integral 
becomes Tm+n = Ge[(m + n)A]T, with: 



'^^rf\Ex\'^ . , „ . 
70 -I sm(26'o) 







(1) 



K being the NLC elastic constant (single constant 
approximation |^) and Aenp the (low frequency) 
anisotropy. When an e-polarized optical beam of slowly- 
varying envelope A propagates in the medium, Eq. 
modifies into: 



KVi 



AeRF\E^\ 



sm{2e) = (2) 



The overall director distribution can be written as 
0{x,y) = do{x,y){l + ip{x,y)), with eo{x,y)ip{x,y) the 
nonlinear optical contribution. In usual experiments Oq 
is small (< 0.4), hence first order approximations are jus- 
tified uf\. The non locality (linked to the V operator in 
Eqs. tias a different impact along x and y, re- 

spectively, due to the strong asymmetry of the problem 
(Fig. 1). A bell-shaped beam with x- waist comparable 
with the cell thickness d does not experience a non lo- 
cal response along x, owing to the planar anchoring with 
^0 ~ in X — zLd/2. Conversely, along y no anchoring 
is present and the index perturbation is free to widen. 
After substituting 9 into Eqs. (E3)-@7 assuming the re- 
sponse to be weakly nonlinear (^ ^ 1) and local in x 
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\fQ{x,y)\ n2^^{x,y)G{C - x) 
xGpirj~y)\fQiC,v)\'dxdydCdrj (6) 



To proceed with the analysis, we need to calculate 
the Green function from Eqs. ^ and J^J). By set- 
ting ix,y) = iXd,YA), 9o = M (^, Y) , ,?(0,0) = 

1, = TM^Hx,n A- - ^^gpp?a^(x,y), 

= // ^(cIx',^Iy) l°(C7??)PdCd?7 (scaling g to § pre- 
vents first order derivatives to appear in Eq. ©); 



R2 _ 'iK a 



^ and 



7(e) = {a + e^^)z?^ + -^ig, equation © with |TJ can 
be cast in the dimensionless form: 

^+/^^ + (i^« + ^/)^ = (7) 



QY2 



7(6)51 + 2^/a^to = 



(8) 



In order to solve Eq. O above, we separate the variables 



3a 
402 



X and Y by letting 
form of Hill's equation in Y: 

dH 



Vy and obtain a simple 



dY^ 



+ {Vy + tf)d = 



(9) 



3 



with Y) = i9{X)^{Y) and i}{X) cx sin(7rX + f ) cor- 
responding to a harmonic oscillator across X. According 
to Floquet theory, the periodic solutions i}{Y) = 'd{Y + l) 
are located on a transition curve !23| . We adopt the per- 
turbative method of strained parameters |2J, [23, per- 
forming the following expansion: 



m 



(10) 



with "dmiY) = ^m{Y + !)■ By collecting and equating 
terms of the same order in e, we find: 



ddi(Yf 

dd2{Yf 
dY^ 



= 

= -Vyld^ - f'&O 



0(1) 

0(6) 



(11) 
(12) 

(13) 



By expanding f{Y) in a Fourier basis / = 
X^m 6xp(— i27rmy), wc obtain the solution to Eqs. 



79=l + ^e",9,„(y) = l+ 

e y exp(-z27rmy) + O(e^) 



i?i(r)/(y)dY + 0(e3 



(14) 



(15) 



Substituting Eq. ifTH) into Eq. (0, for t?(X) « constant 
within the effective modal width, we get: 



[a + ^e™/i„(y)].g + 2V^^/o-0 (16) 



with giX,Y) = uoiX)g{Y), h, = f{Y)[^ - wk)^ + 
2az?i(y) and a periodic /im(F) = h„i{Y + 1) the expres- 
sion of which stems from Eqs. pi|) - (|14|l . For mq = Eq. 

is Hill's equation; otherwise for uq ^ its solutions 
can be found through the expansion: 



5(r) = ^6"g™(r) 



(17) 



Substituting and equating terms of the same order in e, 
we have: 



- ago + 2y/auo 



dg^ 



— -agi^ hi 



90 



^92 h 



0(1) (18) 
0(6) (19) 
0(6^) (20) 



At each order the solution can be factorized in the form 
g{Y) = ge{Y)gpn{Y), i.e. an envelope ge{Y) modulated 
by the periodic function gpniX) — gpn(Y -1-1), with: 



9p0 = 1 

hir 



4:m{a + m^TT^) 



3m exp(— i27rTOK) 



2sign(Y) 



exp(— iTrmy) sin(TO7rK) 



(21) 



Therefore, the generic solution g{Y) can be written in 
the form gniY) = 9eiy)9piY), with a peaked envelope 
gp{Y) and a periodic modulation gp{Y) = ^n^^gpn- 
Now, after introducing the (dimensionless) fields q„ = 

Q^^JZsS^^lMl exp{—i2£^) with ^ — Cz and the degree 
of (non) locality k = yja, by substituting Eqs. (PT|l in 
Eq. Q we finally obtain the discrete non local nonlinear 



Schrodinger equation: 

+ <7n-(-l + <7n-l - 2(7„) 



cothf l^"' 



Equation (|22|l is a discrete version of the non local nonlin- 
ear Schrodinger equation (NNLS) and reduces to the 
discrete nonlinear Schrodinger equation (DNLS) [2^ as 
K — > 00. At variance with previous work dealing with non 
locality in dispersion (a linear property of the medium) 
[l7l I19I I20I I . Eqs. H22|l addresses a nonlinear feature (non 
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FIG. 2: (a) Power Wo of discrete soliton versus size fi for 
increasing k (decreasing non locality) up to the local Kerr- 
case (thick line), (b)-(d) Phase planes of soliton chirp b and 
width jj, for various k in the case jj, = 1.5. 



locality in the nonlinear response) and, thereby, it is ex- 
pected to possess a radically distinct dynamics with re- 
spect to its linear counterpart (see Sec. IIII|) and to non 
local models in the frame of nonlinear photonic crystals 
( [T^ and references therein). The lack of propagation 
terms in Eqs. (8)-(9) of fo^' instance, does not allow 
to study the system evolution (as investigated hereby). 
Hence, equation 1)22(1 can be regarded as a novel general 
model of discrete, dispersive, nonlinear non local media. 



III. ANALYSIS OF THE DISCRETE MODEL 

Equation (|22|) can be also derived using the variational 
principle and the following Lagrangian L: 



I ,2 , exp(-K|m-|-n|) 2, ,2 .oqa 
\qn+i-qn\ +2^ '^^^ih^ 

m 2 



We adopt the peaked ansatz qn{0 = 6xp[«<i5(C) -I- 

ib{^)\n — hqI — fJ-{^)\n — uqW with hq — 0, obtaining the 
effective Lagrangian L^f / : 



dip 



coth /i 



1 db 1 



2 coth fjL 



2 a^ sinh /i 



cosfe .9 , 



(24) 



with non local contribution N{^^ k): 



cxp(— K|m + ?T.|) 



At 



coth(K/2) ' 
tanh(K/2)(2 sinh 2/i + sinhre + sinh(4/i -I- k)) 



(25) 



4 sinh 2/i sin^ (/X + k/2) 
as calculated with a bilateral Z-transform. Setting to 
zero each variational derivative of soliton parameters 
{Aq, </?, b, fi), we obtain the evolution of both soliton chirp 
b and soliton width /x. After some algebra: 



db 4 cos b sinh^ /i 
cosh 2/i 



sinh^ ^ [4iV(/i, k) tanh + 2 sinh^ A^^^^s^] 
cosh 2/x 
4 cosh fi sin & sinh^ fj, 
cosh(2^) 



(26) 



(27) 



with the soliton power Wq — X]„|9n| 



coth /z. 



Equations ((26(1 - ((27|l define a two-dimensional phase- 
space with the conserved Hamiltonian Heff'. 



H, 



ff 



-Wo 



2 + 2 cos 5sech^ + Wq 



tanh(K/2) tanh^ /i(2 sinh 2/i + sinh n + sinh(4/i + k)) 
4 sinh 2/i sin^(/i -|- k/2) 



(28) 



Solitons correspond to stationary points of Eqs. 126(l - l(27() 
with 6 = and: 



Wo 



8 cosh^ fj, cosh I 



3 sinh 



2 



sinh /i sinh^ (/i + ^) 



[sinh/i + sinh(3/i + k)(— 1 -|- 2cosh2/L()] 



(29) 



Since 



tons are stable 



> 0, all fixed points representing soli- 
[23 As shown in Fig. Etj the existence 
curve 1(29(1 of discrete solitons rapidly approaches the lo- 
cal Kerr-case for diminishing non locality (i.e., increasing 
K > 3). As non locality is enhanced and k reduces to- 
wards and below 1, however, the refractive perturbation 
becomes broader and broader (in y) and Wq larger and 
larger. Substantial changes are visible near the soliton 



5 




FIG. 3: (a) Discrete soliton breaking in a quasi-local regime 
(k = 3.0) and (b) breather-like propagation in the non local 
regime (k — 1.0) for (/i — 1.5) and an imprinted chirp 6 = 3. 
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FIG. 4: Calculated reorientation distribution along y {x — 0) 
(solid line) and theoretical Green-function envelope (dashed 
line) in (a) the local (cr = 3.85 x lO^Y/V^) and in (b) the 
non local (cr = 9 X lO^Y/V^) cases, respectively; (c) corre- 
sponding degree of non locality k versus a after a fit (dots) 
and from theory (dashed line); (d) reorientation profile ver- 
sus X {y — 0) without (dashed lines) and with (solid lines) 
optical excitation (1 mW) for local (thick lines) and non local 
responses (thin lines), respectively. 



solution, as displayed in the phase-plane of Eqs. H26() - 
(|27|l in Fig. [2Ib)-(d). In a Kerr regime (k oo) the 
phase-plane consists of a series of periodic orbits near the 
localized state (/i — 1.5, 6 = 0) and fi tends to zero for 
higher chirps b (Fig. [2h)- Therefore, the addition of an 
initial chirp above a certain value -i.e., enough chirp im- 
printing- destroys the soliton 0, IIJ, 28] . The situation 
keeps unchanged as k > 3 (Fig. [5J:). In the non local 
regime {k = 1.0), conversely, the trajectories evolve from 
a closed-loop to a limit-cycle, hence no chirp-imprinting 
can break the soliton (Fig|2i). This remarkable find- 
ing is confirmed by numerical simulations, as visible in 
Fig. 13 While a local system cannot sustain discrete light 
localization with an input spatial chirp above a thresh- 
old (Fig. 13^), non locality allows for the propagation of 
chirped discrete solitary waves (Fig. ^) with periodi- 
cally varying width, as predicted by our model. Clearly, 
the soliton amplitude oscillates as well in order to con- 
serve the total power Wq. These novel solutions belong 
to the class of discrete breathers and are the lattice coun- 
terparts of the continuously breathing solitons reported 
in highly non local bulk NLC 



IV. NUMERICAL RESULTS WITH THE NLC 
LATTICE 

In order to link the analysis to an actual NLC lattice, 
we need to estimate the range of available k. As it stems 
from the model above, the nonlinear index change has a 
peaked envelope of transverse size k = cx AV/^/Kd 
(being Ej.{x) ~ V/d and 9^ ~ const) and, therefore, 
non locality can be tuned by acting on either one of 
the form- factor A/d, the bias V, the elastic constant 
the temperature l^^. With reference to a standard 
NLC (nematic 5CB, with K = i.% y. IO^^^n), to 
evaluate the Green function along Y we employed an 
optical (A = 1.064/.tm) excitation A{x,y) consisting of 
a Dirac distribution across y and a Gaussian beam of 
waist « d across x. We numerically integrated Eqs. 

and Q u sing the relevant potential distribution 
(see Ref. Xh), and finally derived the size n versus 
a = KV/\/Kd by fitting the calculated reorientation 
profile. As the material is tuned from local (Fig. 0^ 
with d/A = 0.44, V = 0.78V) to non local (Fig. Hd with 
d/K — 1.0, y = 0.75V), the Green-function envelope of 
the nonlinear reorientation along y (solid line) is well 
approximated by our theoretical model (dashed line). 
In the latter analysis we employ both geometric {d/K) 
and material {V) tuning, slightly adjusting the bias to 
keep Or ~ 0.35. This results in a quasi-linear transition 
(see Fig. 0]:) from a local to a non local response as 
a varies. The corresponding transverse size of the non 
local response, represented by dots in Fig. shows 
that theory and numerics are in excellent agreement in 
the range covering local (k = 3.0) to non local [n = 1.0) 
responses. Moreover, the all-optical reorientation across 
X, vi sible in Fig. Qlij is nearly sinusoidal (dashed line) 
and does not widen significantly when the nonlinearity 
intervenes (solid line), supporting the validity of the 
local approximation previously adopted. 



V. CONCLUSIONS 

In conclusion, for the first time to the best of our 
knowledge and with specific reference to nematic liq- 
uid crystals and their reorientational response, we have 
modeled discrete light localization in a nonlinear medium 
with an arbitrary degree of transverse non locality. Start- 
ing from the governing equation of the liquid crystalline 
system, we performed an original reduction to a novel 
general form of discrete nonlinear non local Schrodinger 
equation. Remarkably, the latter result was not achieved 
by introducing an a priori specific from of non locality, 
[17, 18, 19, 20J but one derived from the molecular re- 
sponse of NLC. We employed a variational procedure and 
investigated the role of non locality in supporting chirp- 
imprinted discrete spatial solitons. Such novel solutions 
are periodic breathers and cannot exist in purely local 
systems. Since the degree of non locality in NLC arrays 



6 



can be adjusted by acting on geometric or material or 
external parameters p9l |. we anticipate that our findings 
will trigger the observation of discrete light propagation 
in both local and non local regimes in one and the same 
system. Our numerical experiments, in excellent agree- 



ment with the theoretical predictions, fully support such 
possibility. 

We acknowledge enlightening discussions with C. Conti 
and D. Levi. 
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